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The status of the Monte Carlo programs for the simulation of r-lepton production and decay in high-energy 
accelerator experiments is reviewed. No significant changes in the organization of the programs were introduced 
since previous TAU conference, that is why we will concentrate on some physical topics: (i) For TAUOLA Monte 
Carlo generator of r-lepton decays, simulation of five scalar final states based on the hadronic current became 
available for the first time. As an example, simple, but realistic current for final states: 27r~ tv + 2tv° v t , -k~A-k°v t 
and 3ir~ '2ir + v T is presented. The current is installed into TAUOLA. (ii) For the PHOTOS Monte Carlo, which generates 
radiative corrections in arbitrary decays, new results on next to leading order corrections became available for 
some decay modes. The complete corrections were installed for leptonic Z and B decays into a pair of scalars. (iii) 
During conference discussions, the importance of checking the uncertainty of the overall normalization for KORALB 
and KKMC simulations was underlined. Necessary steps to check the uncertainty and to adjust the programs to 
Belle and BaBar conditions are also listed. 
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1. Introduction 

The TAUOLA package |1I2I3I8| for the simulation 
of r-lepton decays and PHOTOS |4I5| for the simu- 
lation of radiative corrections in decays, are com- 
puting projects with a rather long history. Writ- 
ten and maintained by well-defined authors, they 
nonetheless migrated into a wide range of appli- 
cations where they became ingredients of com- 
plicated simulation chains. As a consequence, a 
large number of different versions are presently 
in use. From the algorithmic point of view, they 
often differ only in a few small details, but in- 
corporate many specific results from distinct r- 
lepton measurements. Such versions were mainly 
maintained (and will remain so) by the experi- 
ments taking precision data on r leptons. On 
the other hand, many new applications were de- 
veloped recently, often requiring a program in- 
terface to other packages (e.g. generating events 
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for LHC, LC, Belle or BaBar physics processes). 
The programs organization, prepared for the con- 
venience of users, was presented during previous 
r conference, in year 2004 [6] and we will not re- 
peat it here. 

This time, let us concentrate on more physics 
oriented results. Progress in the simulation of r 
decays into final states of five scalars was achieved 
[7] recently. This work has some consequences for 
the general form how the Monte Carlo programs 
for decay chains have to be organized. We will 
devote section 2 to discussion of that subject. In 
section 3 we will present some new results for the 
simulation of radiative corrections in decays with 
PHOTOS Monte Carlo. Here papers |9|10j . where 
Next to Leading Order effects were introduced 
into generation for the first time, will be pre- 
sented. Even though those results do not affect 
simulation of r decays in the direct way, nonethe- 
less open the way for future precision statements 
on radiative effects in r decays. That is why they 
should find their place in r conference proceed- 
ings. 
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During my talk and in discussions later, I found 
significant interest in domain of overall normal- 
ization of r-lepton pair production cross sections 
as calculated by KORALB [12lllj and KKMC [17 gen- 
erators. These programs were written for higher 
energies. For the use at Belle, BaBar energies 
some effort necessary to adjust photon vacuum 
polarization need to be completed. Section 3 is 
devoted to that purpose. List of necessary bench- 
marks is given and necessary steps are explained. 

Because of the limited space of the contribu- 
tion, and sizable amount of other physically inter- 
esting results, some of them will be excluded from 
conference proceedings. In particular, this time we 
will skip completely the applications related to 
use of r physics in domain of searches for new 
physics. 

2. Five pion final states in TAUOLA Monte 
Carlo 

Approximations were used for five scalar decay 
modes in TAUOLA until recently. In relative recent 
preprint [7], the structure for implementation of 
hadronic currents into r decay modes of five 
scalars was prepared. For the new code, technical 
tests of the solution including benchmark distri- 
bution were performed and described. Finally, 
currents were prepared for some decay modes. 
For the 27r~7r + 27r° mode the three decay chains 
t~ — > a^v — > p~ (— > 7r~7r°)a;(^ 7r~7r + 7r°)i/, 
t~ — > a^v — > ai{-^ 27r~7r + )/o(^ 27r )^, and 
t~ — > a^v — > (— » 7T~27r°)/o(^ 7r + 7r— )v are 
introduced with simple assumptions about the 
couplings and propagators of the various res- 
onances. Similar amplitudes (without the pu> 
contributions) are adopted for the 7r _ 47r° and 
3it~2tt + modes. 

The five-pion amplitude is thus based on a sim- 
ple model, which, however, can be considered as a 
first realistic example. As usual, hadronic current 
is easy to replace by the more sophisticated one, 
once it is required. Also, multitude of additional 
decay channels of five scalars (pions or K-ons) are 
pre-installed. Appropriate flavours and hadronic 
currents have to be provided by the user. 

From that perspective, and in general from the 
point of view of any future Monte Carlo program, 
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27r-7T+27r 


B s. 
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27r-7r+27T° 


A+B s. 


31.35 


31 ±2 


6 


7T~47r° 


B s. 


9.37 


5.5 2 g 


7 


37T-27T+ 


B s. 


11.03 


4.6 ± 0.3 



Table 1: Test results of the generator for real- 
istic choices of parameters; see the text of the 
paper [1] for details. Two currents A and B are 
used in different combinations. The s. (or no s.) 
comment, denotes that for the particular chan- 
nel symmetrization over identical pions was (was 
not) included. With the • we denote the result 
difficult to interpret, because of m^o ^ m^± see 
the text. 



it is of some interest to show some of the numeri- 
cal results. We will concentrate mainly on effects 
of the different types of interferences. 

In the numerical results, collected in Table 
1 the following two currents are used. Cur- 
rent A: r~ — * a^v — ► p~{-^ 7r~7r°)[x)(^ 
7r _ 7r + 7r )^, Current B: r _ — > a^v — > a^(^ 
7T-,7r-7r + (27r ))/ Q (-> 27r (7r + O) i/ - If we com- 
pare lines 3, 4 and 5 of the Table, we can see 
that the interference between current A and B 
seems to be zero (or very small) on total rate. 
It is nice and encouraging result. This picture is 
different when one compare lines 1 with 3, or 2 
with 4. There, interference which appears due to 
symmetrization over identical pions in final state 
is respectively +5 % and -11 %. In case of current 
B, data from the Table are not easy to interpret, 
because of statistical factor | and m^o ^ m^± . 
What is more, the symmetrization effect of -11 % 
is obscured by m T o ^ rn^ effect. 

Let us turn now to more general discussion of 
our numerical results which can be compared with 
the experimental data of ref. [TS] . Perfect agree- 
ment in case of lines 3, 4, and 5 is a trivial con- 
sequence of the fact that these data points were 
used as our input for the current parameters. On 
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the other hand, modest agreement of TAUOLA with 
the data (lines 6 and 7) , provide test of the pre- 
dictive power of our model. 

Rather unpredictive nature of interference ef- 
fects is of importance in construction of Monte 
Carlo programs, if they are supposed to be used 
in comparisons with the precision data on decays. 
The constructive/destructive interferences origi- 
nate from the corrections of order T/Q for in- 
termediate resonances, which appear in the cas- 
cade decays. If widths of the resonances satu- 
rating currents would be sufficiently small then 
of course interferences would be negligibly small. 
At the same time the intermediate states would 
be well formed and necessities to use different 
parametrizations for example of p resonance, de- 
pending on how it was formed and how it will 
decay, would not be necessary. 

Unfortunately, in practice, it is not the case, 
and we have to bear this constraint unpleasant 
for program construction in mind. 

3. PHOTOS and NLO effects in B and Z de- 
cays 

There were significant changes introduced into 
PHOTOS Monte Carlo project over the last two 
years. The complexity of the subject matches 
neither size nor the purpose of the present talk. 
Recently collected results [9110] only indirectly af- 
fect predictions of simulation for r physics. How- 
ever as for the first time explicit form of the 
approximate matrix elements actually used was 
written, new possibilities opened. In particu- 
lar, for the two-body decay modes into fermions 
[Z — > or scalars (eg. B — » K + tt~) ex- 

act NLO order matrix elements were implemented 
into PHOTOS kernel. In cases when the correction 
kernels were switched on and simulations were re- 
stricted to the first order, the differences with re- 
sults of reference Monte Carlo programs were be- 
low statistical errors of 10 8 event samples. Also 
technical tests of PHOTOS performed at that or 
even better statistical level confirmed the correct- 
ness of program design; even though significant 
changes in crude distributions were introduced, 
the results did not change at all. 

To visualize the results let us present the ex- 



ample plots for (Z —> and B — > K + tt~ 

cases. In case of (Z — > fi + [i~) we could perform 
comparisons for multi-photon version of the gen- 
eration, because our reference program KKMC has 
such possibility as well. The histogram with the 
worst agreement of all possible invariant masses 
which could be constructed from four momenta 
of: [i~ and eventually up to two hardest pho- 
tons (their energies had to be above 1 GeV) is 
shown in Fig 1 where the invariant mass of 77 
pair is plotted. Even in this case, the green and 
red (gray) lines for PHOTOS and KKMC nearly over- 
lap (the logarithmic left side scale has to be used 
for that lines). To visualize the differences, we 
show on the plot the ratio of the two histograms 
which is depicted with the black line and follows 
linear right scale. The differences in the ratio are 
rather small for the masses of 77 pair up to about 
50 GeV. For larger masses which contribute about 
10 -3 to the sample of two-photon configurations 
the difference gradually grow to about 15% close 
to the phase space limit. Note that fraction of 
events with least two photons each of which with 
energies above 1 GeV is 1.26 % from KKMC and 
1.29 % from PHOTOS. 

In case of the B meson decays the agreement 
between PHOTOS and reference calculation, even 
without the use of correcting weight was excel- 
lent, that is better than for the presented above 
result for Z decay. That is why we skip numer- 
ical result and address the reader to conference 
transparencies or to ref. [9]. 

4. Normalization issues for Monte Carlo 
programs for r-lepton pair production 
at Belle/BaBar energies 

In discussions it became obvious that some 
comments on the overall normalization of KORALB 
and KKMC predictions for r-lepton pair production 
at Belle and BaBar energies is of importance for 
the present day users. 

One has to bear in mind that KORALB was pub- 
lished [12111] more than twenty years ago. The 
program was supposed to feature orthodox first 
order QED corrections and complete mass and 
spin effects. Such formula turned out to be very 
useful, and program remains in broad use until 
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Figure 1. Comparisons of improved PHOTOS with mul- 
tiple photon emission and KKMC with second order ma- 
trix element and exponentiation. Comparison figure 
of worst agreement was selected from all constructed 
with the help of MC-TESTER fig)? for details see the text 
of the paper. 



now. On the other hand, some of its inputs are 
rather outdated and do not match the present 
day requirements, even for technical tests. To be 
precise, I have in mind the function PIRET(S), 
which features the real part of hadronic vacuum 
polarization of photon as measured by the data 
collected until early 80 's. Obviously this function 
need to be replaced by the more modern one. A 
possible choice can be the function REP I of ref. 
[T5] because it t has similar functionality as PIRET 
of KORALB. 

Unfortunately the improvements on FORTRAN 
function PIRET (S) describing hadronic vacuum 
polarization of photon, do not solve all normal- 
ization problems of KORALB. It is well known since 
long, that the genuine one loop corrections are not 
enough, and the solutions are available. The two 
major improvements which appeared as a conse- 
quence of phenomenology improvements during 
the LEP era was introduction of higher order 
QED corrections into Monte Carlo simulations 
and better way to combine loop corrections with 
the rest of the field theory calculations. It was 



found to be safe to sum contributions of loop cor- 
rections into photon (and Z) propagators. Then, 
terms of all, but incomplete, orders of perturba- 
tion expansion are taken into account. That is 
why significant effort was needed to justify the ap- 
proach [14] . At lower energies things are of course 
simpler, as there is no need to worry about Z- 
lineshape. The KKMC Monte Carlo [17] could thus 
be a complete solution to B-factory needs and 
ready to use. This program features higher order 
QED matrix elements with the help of exclusive 
exponentiation, and explicit matrix element up 
to the second order. Unfortunately electroweak 
library (which include vacuum polarization, the 
only function interesting for low energy applica- 
tions) was never adopted to requirements of pro- 
cesses below 13 GeV center-of-mass energies. 

To provide reasonably good results, this library 
needs thus to be overruled (analogously to the 
function PIRET (S) of KORALB) by a more suitable 
oncQ. After small tests the function REPI of ref. 
13J should be suitable (overall normalization con- 
stants and other details of conventions need to be 
checked). 

Once this is completed, and if two-loop photon 
vacuum polarization can be neglected, KORALB 
and KKMC can form a base for tests and studies of 
systematic errors for cross section normalizations 
at low energies. The necessary strategy should be 
similar to the one for the Bhabha scattering see 
eg. [TH] for description. 

For that purpose the following calculations of 
the total cross section should be performed: 

1. KORALB radiative corrections switched off. 

2. KKMC, both electroweak corrections and ini- 
tial state QED bremsstrahlung switched off 
(final state bremsstrahlung need to be kept 
on, otherwise spin amplitudes necessary for 



1 In KKMC similar improvements require to overwrite in 
routine DZf ace_MakeGSW (file dizet/Dzf ace . f ) calculation 
of GSW(6) with the one using function PIRET(SS) or 
REPI(SS). The PIRET(SS) must be supplemented with the 
leptonic contribution to vaccum polarization as well. It 
was checked by S. Banerjee that pretabulation (in contrast 
of low energy numerical values of KKMC vaccum polar- 
ization) is not a problem. To this end he simply increased 
density of pretabulation points. 
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calculation of r spin correlations are not cal- 
culated and program stops). 

3. Radiative corrections switched on in 
KORALB. 

4. Radiative corrections switched on in 
KORALB. Vacuum polarization switched off 
with the help of internal key IFVAC=0. 

5. KKMC: electroweak corrections on, initial 
state QED bremsstrahlung switched off. 

6. KKMC: electroweak corrections switched off 
but initial state QED bremsstrahlung 
switched on 

7. KKMC: both electroweak corrections and ini- 
tial state QED bremsstrahlung switched on 

The results, we will call them respectively 
cr±... (77, can be calculated for r- or /i-pair pro- 
duction. Also experimental cuts can be applied 
in the calculation of these cross section. Let us 
list now some of the possible checks. Of course all 
calculations being compared must be performed 
with the same assumptions on experimental cuts 
and final state flavours. 

• : g 2 should be larger by a factor of 1 + |- 
than <7i for all center of mass energies in case no 
experimental cuts are applied. •: the following 

relation^ should hold: a 3 - cr 4 = 0"i2(l - y^ff)- 
• : finally comparisons of cr 6 a 4, with a\ can be 
used to estimate the size of QED bremsstrahlung 
effects respectively of higher orders or just first 
order alone. •: the relation — = ^ is not ex- 
pected to hold precisely. It can hint on numerical 
importance of convoluting QED bremsstrahlung 
and vacuum polarization with respect to naive 
factorization. 

If comparisons are repeated with experimen- 
tal cuts applied, some extra care must be taken, 

2 In case of KORALB vacuum polarization increase the cross 
section by (25ftll 7i7 (s)) times Born cross section. In case of 
KKMC the vacuum polarization factor 1 rpjp appears 

for all events. In case when initial state bremsstrahlung 
is switched on, the effective transfer s' depends on the 
amount of energy carried out by the initial state radia- 
tion convolution with initial state radiation spectrum is 
formed. 



because of cut off dependence of final state 
bremsstrahlung effects. On the other hand, nu- 
merically significant and theoretically ambiguous 
contributions from events with very low final state 
lepton pair invariant masses, are removed. 

Unfortunately the above points must remain in 
sketchy form. Full clarification, as LEP experi- 
ence showed, require significant amount of work. 

5. Summary and future possibilities 

The status of the computer programs for the 
decay of r leptons and associated projects was re- 
viewed. The high-precision version of PHOTOS for 
radiative corrections was presented. In particu- 
lar, the option to run the program with multiple- 
photon radiation was mentioned. New results for 
leptonic decays of Z, and B meson decays into 
pair of scalars, were presented. For these channels 
complete next-to-leading order effects can now be 
simulated. However, for the most of the applica- 
tions these effects are not necessary, leaving stan- 
dard modular version of PHOTOS sufficient. The 
important result of the above work, is that the 
path to include electromagnetic form-factors of 
the particles participating in decay is now open 
for the future fits to the data. These form-factor 
effects may be significantly larger and physically 
more justified, than complete next-to-leading or- 
der effects of scalar QED in B meson decays re- 
cently installed,, 

The presentation of the TAUOLA general- 
purpose interface was omitted because of lack of 
time. Examples for its use in the case of the Higgs 
boson parity measurement at a future Linear Col- 
lider [19120121] and for Higgs searches at the LHC 
[22"] can be found in the literature. Recently, a 
similar application was developed for the case of 
studies in hypothetical effects of CP-parity break- 
ing in the Bq-Bq system at Belle and BaBar [23 . 

Distinct versions of the TAUOLA library for r 
lepton decay, and of PHOTOS for radiative correc- 
tions in decays, are now in use. The principles 
how to use the distribution package, are presented 
inrefs. [SIS] , 

In the talk we have reviewed the results for the 
novel decay modes of TAUOLA into five scalars. 
These modes feature now the hadronic current. 
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The simple but realistic current [7] is avail- 
able for decay modes 27r _ 7r + 27r°;/ T , tt^Att ^ and 
3it~2tt + v t . Numerical study of the new decay 
modes helped to formulate comments on the im- 
portance of T/M terms for the intermediate res- 
onances. It is argued that the constructive or de- 
structive interferences appear and that their ex- 
istence must be taken into account by builders of 
the future Monte Carlo programs for decays (not 
necessarily r decays). 

Finally, presentation of adjustments for KORALB 
and KKMC programs in treatment of photonic vac- 
uum polarization was given. It was explained, 
that without such changes the programs can not 
be used for discussion of the normalization un- 
certainty for the r-lepton pair production cross 
section at Belle/BaBar energies. 

Acknowledgements 

Useful collaboration and suggestions from E. 
Barbcrio, G. Nanava, S. Jadach, J. H. Kiihn, 
S. Eidclman, B. Kersevan, E. Richter-Was P. 
Golonka are acknowledged. Discussions with 
members of the Belle and BaBar collaborations 
are also acknowledged. Finally, exchange of e- 
mails with B. Pietrzyk, H. Burkhardt, D. Bardin, 
T. Riemann, S. Banerjee and M. Roney was a nec- 
essary input to section 4 of the present report. 

REFERENCES 

1. S. Jadach, J. H. Kuhn, and Z. Was, Comput. 
Phys. Commun. 64 (1990) 275. 

2. M. Jezabek, Z. Was, S. Jadach, and J. H. 
Kuhn, Comput. Phys. Commun. 70 (1992) 
69. 

3. S. Jadach, Z. Was, R. Decker, and J. H. Kuhn, 
Comput. Phys. Commun. 76 (1993) 361-380. 

4. E. Barbcrio, B. van Eijk, and Z. Was, Com- 
put. Phys. Commun. 66 (1991) 115. 

5. E. Barberio and Z. Was, Comput. Phys. Com- 
mun. 79 (1994) 291-308. 

6. Z. Was and P. Golonka, Nucl. 
Phys. Proc. Suppl. 144 (2005) 88 

arXiv:hep-ph/0411377l . 

7. J . H. Kuhn and Z. Was, 
arXiv:hep-ph/0602 1621 

8. P. Golonka et at, Comput. Phys. Commun. 



174 (2006) 818 [arXiv:hep-ph/0312240] . 

9. G. Nanava and Z. Was, 
|arXiv:hep-ph/06070"l9 

10. P. Golonka and Z. Was, 
|arXiv:hep -ph/0604232] 

11. S. Jadach and Z. Was, "QED O (Alpha**3) 
Radiative Corrections To The Reaction E+ 
E- -> Tau+ Acta Phys. Polon. B 15 (1984) 
1151 [Erratum-ibid. B 16 (1985) 483]. 

12. S. Jadach and Z. Was, "Monte Carlo Simu- 
lation Of The Process E+ E- -> Tau+ Tau- 
Including Comput. Phys. Commun. 36 (1985) 
191. 

13. H. Burkhardt and B. Pietrzyk, Phys. Rev. 
D 72 (2005) 057501 [arXiv:hep-ph/0506323| . 
Fortran routine REPI(S), 
|http: / /hbu.home.cerrL ch/hbu/aqed /repi2005 

14. D. Y. Bardin and G. Passarino, "The stan- 
dard model in the making: Precision study of 
the electroweak 

15. S. Eidelman et al. [Particle Data Group Col- 
laboration], Phys. Lett. B 592 (2004) 1. 

16. P. Golonka, T. Picrzchala, and Z. Was, 
Comput. Phys. Commun. 157 (2004) 39-62, 
|hep-ph/0210252| 

17. S. Jadach, B. F. L. Ward and Z. Was, "The 
precision Monte Carlo event generator KK for 
two-fermion final states Comput. Phys. Com- 
mun. 130 (2000) 260 |arXiv:hep-ph/9912214| . 

18. S. Jadach, E. Richter-Was, B. F. L. Ward and 
Z. Was, Phys. Lett. B 260 (1991) 438. 

19. G. R. Bower, T. Pierzchala, Z. Was, and 
M. Worek, Phys. Lett. B543 (2002) 227-234, 
|hep-ph/0204292[ 

20. K. Desch, Z. Was, and M. Worek, Eur. Phys. 
J. C29 (2003) 491-496, |hep-ph/0302046l 

21. K. Desch, A. Imhof, Z. Was, and 
M. Worek, Phys. Lett. B579 (2004) 157-164, 
|hep-ph/0307331| 

22. E. Richter-Was, T. Szymocha, and 
Z. Was, Phys. Lett. B589 (2004) 125- 
134, hep-ph/0402159| 

23. P. H. Chankowski, J. Kalinowski, Z. Was and 
M. Worek, Nucl. Ph ys. B 713 (2005) 555 
[arXiv:hep-ph/0412253] . 



